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GRAVITY WATER FLOWS WITH DISCONTINUOUS VORTICITY AND 

STAGNATION POINTS 

CALIN IULIAN MARTIN AND BOGDAN-VASILE MATIOC 


Abstract. We construct small-amplitude steady periodic gravity water waves arising as 
the free surface of water flows that contain stagnation points and possess a discontinuous 
distribution of vorticity in the sense that the flows consist of two layers of constant but 
different vorticities. We also describe the streamline pattern in the moving frame for the 
constructed flows. 


1. Introduction 

We present here a study of steady periodic traveling water waves that propagate at the 
free surface of a two-dimensional inviscid and incompressible fluid of finite depth allowing 
for stagnation points and for a discontinuous distribution of vorticity. More precisely, we 
consider water waves interacting with two vertically superposed currents of different constant 
vorticities. 

Confined first to the investigation of waves of small amplitude, which can be satisfactorily 
approximated by sinusoidal curves within the linear theory, the examination of periodic 
traveling water waves arising as the free surface of an irrotational flow with a flat bed 
originates at the beginning of the \Q th century. The description of waves that are flatter 
near the trough and have steeper elevations near the crest necessitates a nonlinear approach, 
which was in fact conducted in the last decades and led to the first rigorous results concerning 
the existence of wave trains in irrotational flow, see for instance the case of Stokes waves [34] 
and the flow beneath them (particle trajectories, behavior of the pressure) cf. [2, 3, 6, 8]. 

To go beyond irrotational flows and to treat wave current interactions one needs to in¬ 
corporate vorticity into the problem, cf. [4, 16, 32]. However, the difficulties generated by 
the presence of the vorticity have prevented a rigorous mathematical development, which 
appeared only relatively recently in [7], where the existence of small and large amplitude 
steady periodic gravity water waves with a general (continuous) vorticity distribution was 
proved. 

Of high significance is the investigation of steady periodic rotational waves interacting 
with currents that possess a rough - that is discontinuous or unbounded - vorticity. Discon¬ 
tinuous vorticities model sudden changes in the underlying current, numerical simulations of 
such flows being quite recent [18, 19]. Unbounded vorticities on the other hand can describe 
turbulent flows in channels (see the empirical law in on page 106 of [1]) and are relevant also 
in the setting of wind generated waves that possess a thin layer of high vorticity adjacent 
to the wave surface [30, 31]. The discontinuous vorticity distribution was considered in the 
groundbreaking paper [9] where the existence of steady two-dimensional periodic gravity 
water waves of small and large amplitude on water flows with an arbitrary bounded (but 
discontinuous) vorticity was proved. Small amplitude capillary-gravity waves with discon¬ 
tinuous but bounded vorticity were constructed in [25, 29]. Waves with unbounded vorticity 
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were first shown to exist in [26] but only when allowing for surface tension as a restoring 
force. This situation appears in many physical settings one of which being that of wind 
blowing over a still fluid surface and giving rise to two-dimensional small amplitude wave 
trains driven by capillarity [17] which grow larger and turn into capillary-gravity waves. 

Another striking occurrence in water flows is the presence of stagnation points, that is 
points where the steady velocity field vanishes, thus making the analysis more intricate, since 
the usual Dubreil-Jacotin transform which converts the original free boundary problem into a 
problem in a fixed domain, is no longer available. There is a short list of papers dealing with 
existence of water flows allowing for stagnation points and for a non-vanishing continuous 
vorticity, cf. [10, 13, 14, 20, 35] for gravity waves and [22, 23, 24, 28] for waves with 
capillarity. Under consideration in this paper is a a more involved setting where, in addition 
to permitting stagnation points (whose existence in the fluid is proven), we also allow for 
a discontinuous distribution of the vorticity. To our best knowledge the incorporation of 
both stagnation points and of a discontinuous vorticity is a feature that was not rigorously 
analyzed before. 

The governing equations are the Euler equations of motion, together with boundary 
conditions on the free surface and on the flat bed of the water flow. The discontinuous 
vorticity that we consider here is of the following type: we assume that the flow has a layer 
of constant vorticity 72 adjacent to the free surface above another layer of constant vorticity 
71 neighboring the flat bed. Of course, the interesting situation (that we pursue here) is 
when 71 7 ^ 72 - The unknowns are here the free surface, the interface separating the regions 
of different vorticities (which can be seen as an internal wave due to the discontinuity in 
vorticity), the velocity field and the pressure function. In a first step we reduce the number 
of unknowns by means of the stream function whose utilization converts the problem into 
a transmission problem along the line of discontinuity of vorticity with fewer unknowns. 
The second step that we undertake is to consider a flattening transformation which has the 
advantage that changes the free boundary value problem into a problem in a fixed domain, 
thus making it more tractable for the analysis. For studying the latter resulted problem we 
employ the Crandall-Rabinowitz Theorem on bifurcation from simple eigenvalues. 

The dispersion relation that we obtain - which is a formula giving the speed at the 
free surface of the bifurcation inducing laminar flows in terms of the two vorticities 71 , 72 , 
the thickness of the two rotational layers and the wavelength - generalizes the one in [ 10 ] 
obtained in the case of a water flow with constant vorticity and allowing for stagnation 
points. The intricacy of the dispersion relation - a third order algebraic equation - allows 
us to prove existence of water waves of small wavelength arising as the free surface of water 
flows with rotational layers of different constant vorticities and containing stagnation points, 
cf. Theorems 3.4-3. 6 , 4.2, 4.3. We present also the streamline pattern in the moving frame 
for the solutions, cf. Figures 1-3. Our results show especially that the ratio of the amplitudes 
of the surface wave and that of the internal wave - and the fact that the surface wave and 
the internal wave are in phase or anti phase - is highly influenced by the vorticities of the 
currents and by the speed at the free surface of the bifurcation inducing laminar flows. 

We briefly outline the content of the paper. We present in Section 2 the governing equa¬ 
tions together with the analytic setting we work in. Moreover, we also find the dispersion 
relation whose analysis is undertaken in Section 3 for the case 72 > 0, while the more singular 
case 72 = 0 is treated in Section 4. The Appendix contains several technical lemmas. 
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2. The governing equations 


Under consideration is a two-dimensional steady periodic flow, moving under the influence 
of gravity, such that the surface waves propagate in the positive x-direction. The water flow 
occupies the domain U bounded below by the flat bed y = —d, with d > 0, and above by 

the free surface y = h(x), which is a small perturbation of the flat free surface y = 0. In a 

reference frame moving with the wave speed c > 0, the equations of motion in II are Euler’s 
equations 

(u - c)u x + vu y = -P x , 

< (u - c)v x + vv y = -Py-g, (2.1a) 

^ Hx T U/ — 0) 

where (u,v) denotes the velocity field, P stands for pressure and g is the gravity constant. 
The equations of motion are supplemented by the boundary conditions, which, ignoring 
surface tension effects, are 

{ P = P 0 on y = h(x), 

v = (u — c)h' on y = h(x), ( 2 . 1 b) 

v = 0 on y = —d, 

with P 0 being the constant atmospheric pressure. 

We are interested in solutions of the problem (2.1) for which the vorticity u := u y — v x 
of the flow presents discontinuities of the following type: we assume that, adjacent to the 
free surface, the water flow possesses a layer 

^(/, h) := {(x, y) : x € R, -d 2 + f(x) < y < h{x)}, 

of constant vorticity 72 , situated above another layer 

^(/) : = {Ou y) ■ x €R,-d < y < -d 2 + f(x)}, 


which is adjacent to the flat bed and is of constant vorticity 71 , that is 


/ 7i) in ^(/)> 

1 72, in 


( 2 . 1 c) 


We steadily assume that 71 7 ^ 72 and that d 2 > 0, d — d 2 =: d\ > 0. We note that, 
additionally to (u, v, P, h ), we have a further unknown: the function / whose graph separates 
the two currents of different vorticities. By Helmholtz’s law, the vorticity is constant along 
streamlines of the steady flow, and as a consequence of this y = —d 2 + f has to be a 
streamline of the flow. This streamline can be viewed as an internal wave due to the jump 
in vorticity. 

With the help of the stream function ip, introduced (up to an additive constant) via the 
relation Vi/> = (—v,u — c) we can reformulate ( 2 . 1 a)-( 2 . 1 c) as the free-boundary problem 


subjected to the conditions 


A?/>2 = 72 

in H(/, h), 

= 71 

in n(f), 

Ip 2 = 0 

on y = h(x), 

1p2 = Ipl 

on y = —d 2 + f(x), 

ipi = m 

on y = —d, 

■ dylpi 

on y = -d 2 + f(x 


| VV^ 2 1 2 + 2g(d T h) = Q on y = h(x), 


( 2 . 2 a) 


( 2 . 2 b) 
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where the constant — m represents the relative mass flux and Q € R is related with the 
hydraulic head. Moreover, i/q := i/j and il> 2 := V 7 1si(/ft.)’ so that ^ rom the f° ur th equation 
of (2.2a) and the first equation of (2.2b) we see that the \7ip (hence also the velocity field) 
is continuous across the interface y = —c ?2 + f( x )- 

Given a € (0,1), it is not difficult to see that any solution 

e (c p 3 +“(m )) 2 X c 3 + a (n(f)) x c 3 p + a mTh)) 

of ( 2 . 2 ) defines a solution 

(u,v,P, ( f,h )) € {C p ~(Ti )) 3 x (C 3 + a (R )) 2 _ 

((“> T, )ln(f)>(“> t, )lnff,h)) G ( C per a {V(f ))) 2 X (C p 2 +“(H(/,h ))) 2 

( p \n U y p Wm) € C per a (^f)) x 

of (2.1). The subscript per stands for functions that are periodic in the horizontal variable, 
meaning that all the functions considered above are L— periodic with respect to x, with 
L > 0 being fixed. 

We first determine laminar flow solutions of problem (2.2), that is water flows with 
a flat free surface and parallel streamlines, meaning that they present no x-dependence. 
Of interest are laminar flows that contain stagnation points, more precisely laminar flows 
that contain streamlines consisting entirely of stagnation points. Then we study when non 
laminar solutions bifurcate from the laminar flows and describe the qualitative picture of 
the streamline pattern for the bifurcating solutions. 

Laminar flow solutions. Because the stream function is constant along the streamline 
y = — d 2 + f{x), we use the value of the stream function 

V’l = 4>2 = A on y = -d 2 + f(x), (2.3) 


to parametrize a family of laminar solutions of (2.2a). Setting / = h = 0 we obtain from 
( 2 . 2 a) that the stream function i/>° := ( 1 ^ 1 , ^ 2 ) satisfies 
,2 


W(y) = 


7i 2T 


+ 


71 (d + d 2 ) A - m 


2 ' di 

, 0 , ^ 72l/ 2 , fl2d 2 \\ _ r j nl 

Mv) = — y 6 | - <i2 '°]' 

The equations of (2.2b) are equivalent to 


A d 71 dd 2 md 2 

y+ i; + —— ~dT’ 


A d 7 ihi + 72^2 

m = — + di --- ; 

a 2 Z 


Q = 


n_2i2_\y + 2gi 


\ 2 


d 2 J 


(2.4) 


In the following we choose the constants m and Q in (2.2a) and (2.2b) to be given by (2.4), 
the constant A introduced via (2.3) being left as a parameter. Hence, each A £ R determines 
a unique laminar solution ((f,h),ijji,'ip 2 ) := (0,i/^) of (2.2) when m and Q are defined 
by (2.4). These are the laminar solutions from which we study bifurcation. 

Conditions for stagnation. We note that the laminar flows determined above possess 
stagnation points - that is water particles that travel horizontally with the wave speed - if 
and only if 


d y ip 2 (-d 2 ) • dyip® ( 0 ) <0 or d y i />°(—d) • d y il^_{—d 2 ) < 0 . 


If (2.5) holds true, then there exists yo £ [— d, 0] such that 

d y ipi(yo) = 0 for i = 1 or 2 . 


(2.5) 
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The streamline y = yo consists only of stagnation points, and we expect that the solutions of 
(2.1) that bifurcate from these laminar solutions possess stagnation points too, cf. [12, 35]. 
The first inequality ensures stagnation in the layer adjacent to the wave surface, and is 
equivalent to 

A(A - 7 2 d 2 ) < 0 , ( 2 . 6 ) 

respectively the condition for stagnation in the bottom layer is 

(A - 7 2 d 2 )(A — 7 idi - 72 ^ 2 ) < 0. (2.7) 


Hereby, we set 


A : = 


72^2 


A_ 

T 2 


( 2 . 8 ) 


The constant A has a physical interpretation: it is the relative horizontal speed at the free 
surface for the laminar flow determined by A, that is A = dyijj®\ y ~Q- For this reason we 
define A via (2.8) and use A as parameter. 

The analytic setting. With A as parameter, we are left to seek special values of A such 
that branches of non laminar solutions of ( 2 . 2 ) bifurcate from the curve of laminar flows. 
For this, we need to recast (2.2) in a suitable analytic setting. 

In the following a € (0,1) is a fixed Holder exponent. Because the equations of (2.2a) and 
( 2 . 2 b) are posed on manifolds that depend on the unknown functions (/, h), it is suitable to 
transform the problem (2.2) on fixed manifolds. For this, we set 12i := 11(0), 12 2 '■= 12(0,0) 
and define the mappings 


: 17i —> 12(/), $ f (x, y) = (x, dl + d ^ V + ^ '~f (z)), 

$(. f,h) ■ ^2 -t 12 (f,h), $ (M) (x,y) = (x,^^- - /j x ) + d2 y + h( x )J . 

It is easy to see that and are C' 3+ “ —diffeormorphisms for each (/, h ) € O, whereby 


0:={(f,h)€ (Cj+£(R )) 2 : -d < -d 2 + f < h}, 

the subscript e referring to the fact that we consider only even function in x. Using these 
diffeomorphisms, we define the linear elliptic operators 

A(f) : CfgM) Cl+Z(U 1 ), A(f) Wl := A(uq o &J 1 ) o 4> /; 

A(f, h ) : C 3 +“ (fi 2 ) -a 0(12 2 ), A(f, h)w 2 := A {w 2 o ^ h) ) o <h (M) , 
and the boundary operators 

Bi : M x O x Cj+£(n 2 ) -A C 2 + Q r (M), 

B 2 :RxOx C 3 +“ (HO x C 3 +“ (12 2 ) -a C 2 +“ (M), 


respectively through 

£u(A, (/, h),w 2 ) := (\V(w 2 o$-* h) )\ 2 o$( f>h) +2g(d + h) - q) 


y =0 


B 2 (A, (/, h))[w 1 ,w 2 ] 


(dy ( W2 O $ U ] h) )) O $ {m - (dy(w 1 O $/)) O 


Observation 2.1. Let A € M, ((/, / 1 ), ipi, ^ 2 ) € 0 x C 3 + Q ( 12(/)) x C 3 +“( 12(/,/i)), and 
assume that \,m,Q are defined by (2.3), (2.4), and (2.8). Then, the tuple ((/, /i), 1 / 7 , i/> 2 ) 
solves the problem ( 2 . 2 ) if and only if 











6 


C. I. MARTIN AND B.-V. MATIOC 


(*) Wl ■= V’l ° (fii) is the unique solution of the Dirichlet problem 

{ _4.(/)xdi = 71 in fii, 

uq = A on y = — d 2 , (2.9) 

w\ = m on y = —d. 

(ii) w 2 '■= ip 2 ° € C*g p^(n 2 ) is the unique solution of the Dirichlet problem 

{ A(f,h)w 2 = ^ 2 in D 2 , 

w 2 = 0 on y = 0 , ( 2 - 10 ) 

w 2 = A on y = —d 2 . 

(m) £i(A, ( f,h),w 2 ) = # 2 (A, (/,fr))[wi,W 2 ] = 0 in 

Thanks to Observation 2.1 we can recast the problem (2.2) as a nonlinear and nonlocal 
equation with (A, ( f,h )) as unknown. In order to proceed, we establish first the following 
result. 

Lemma 2.2. Given (A, (/, h )) € R x O, we let w\ := wi(A, (/, h )) and w 2 '■= w 2 ( A, (/, /i)) 
denote the unique solution of (2.9) and (2.10), respectively, with A given by (2.8). T/ien, 
we have Wi € C^M x C>,C'|+^.(fij)),i = 1,2. 

Proof. We prove just the real-analyticity of the solution operator w i, the claim for u > 2 
following similarly. By elliptic theory, cf. e.g. [15], we see that w\ : R x O —>• C'gp^-(Di) is 
well-defined. Moreover, we have that 

^(A, (/, h),w i (A, (/, /!))) = 0 for all (A, (/, fc)) g I x O, 

whereby JsC“(lxOx Ce,per(^i), ^e,per(^i) x (^e,perW) 2 ) i s the operator defined by 

T(A, (f,h),w 1 ) := (A(f)wi - 'y 1 ,w 1 \ y= _ d2 ,w 1 \ y= _ d ). 

Taking into account that Frechet derivative 

d Wl iF(A, (/, h),wi((A, (/, h))))[z] := {A(f)z, z\ y= _ d2 , z\ y= _ d ) 

is an isomorphism, the assertion follows from the implicit function theorem. □ 

Because £>,;,* = 1,2, depend real-analytically on their arguments too, we obtain from 
Lemma 2.2 and Observation 2.1 that the problem (2.2) is equivalent to the nonlinear and 
nonlocal equation 

4>(A,(f,h)) = 0, (2.11) 

whereby <h := ($i,<h 2 ) € C u (R X O, (Cgjj“ r (M)) 2 ) is the operator defined by 

<K(A, (/, h)) := (Bi(A, (/, h),w 2 ( A, (/, h))), B 2 { A, (/, /i))MA, (/, fc)), w 2 (A, (/, />))]). 

( 2 . 12 ) 

The laminar flow solutions of ( 2 . 2 ) correspond to the trivial solutions (A, 0) € R x O of 
(2.11). In order to find other solutions, we use the theorem on bifurcations from simple 
eigenvalues due to Crandall and Rabinowitz [11], 

Theorem 2.3 (Crandall and Rabinowitz). Let X, Y be real Banach spaces and let the map¬ 
ping $ € C U (WL x X, Y) satisfy: 

(a) <h(A, 0) = 0 for all A € R; 

(b) There exists A* € R such that Frechet derivative <9 X <1?(A*,0) is a Fredholm operator 
of index zero with a one-dimensional kernel and 

Ker (9a;$(A*, 0) = spanjxo} with 0 7 ^ xq £ X; 
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(c) The transversality condition 

<9a x $(A*,0)[x 0 ] £ Im<9 x 4>(A*,0). 

Then, (A*,0) is a bifurcation point in the sense that there exists £ > 0 and a real-analytic 
curve (A, x) : (—•£,£) -} IxX consisting only of solutions of the equation 4>(A, x) = 0. 
Moreover, as s —>• 0, we have that 

A(s) = A* + O(s) and x(s) = sxo + 0(s 2 ). 

Furthermore, there exists an open set U ClxX with (A*, 0) G U and 

{(A, x) G U : 4>(A,x) = 0,x / 0} = {(A(s),x(s)) : 0 < |s| < e}. 

In order to apply this abstract bifurcation result, we need to compute the Frechet deriv¬ 
ative of the operator 4>. To this end we state the following lemma. 

Lemma 2.4. Let A G R be given. The Frechet derivative ^(/,/i)^(A, 0) is the matrix operator 
d m m,0) = (A; F 2 ) € £((Cj«(K)) 2 , (C?+“(R)) 2 ). 


symbol (A) := (mV (A))^^ defined by 


\OL 

i,per 


Given 1 < i, j < 2, the operator Afj G C (C^ t 
ij 
l k 

11 / A \ a A / . .7 A \ 


/^ 2 +. 

5 


Q! 

e,per 


m k (A) = - 2 A( 7 2 d 2 - A)- 

m fc 2 ( A ) = 2^ + 72A - A 


sinh (Rk/h ) ’ 

2 


ml 1 (A) = 72 - 7l + (A - 72 ^ 2 ) 


tanh (it^efo) J ’ 
Rk 


+ 


-tanh (R k d\) tanh (R k d 2 )i ’ 


(A) = -A- 




zs t/ie Fourier multiplier with 

(2.13) 

(2.14) 

(2.15) 

(2.16) 


Rk 


sinh ( R k d 2 ) 

/or k G N, whereby R k := 2kTr/L. For k = 0 f/ie right-hand side of (2.13)-(2.16) should be 
understood as the limit of the expressions when letting R k —X 0. 

Proof. See Appendix. □ 


With the help of Lemma 2.4 we are now able to determine when the Frechet derivative 
d(f,h) ^(A, 0) is a Fredholm operator. 

Lemma 2.5. Let A G M be given. We have: 

(i) If A G { 0 , 72 ^ 2 }) then A, 0) is not a Fredholm operator. 

(■ ii ) If A 0 { 0 , 72 ^ 2 }) i/ien 0) is a Fredholm operator of index zero. 

Proof. In order to prove (i), we infer from (2.13) and (2.14) that for A = 0 we have 

d {f , h) $i(A, 0)[(/, h)} = 2 gh for all (/, h) G (C 3 +“(M)) 2 , 

meaning that Im9(j^$i(A,0) = C 3 +“ r (M). Since C 3 +^,(M) is not a closed subspace of 
Cgp^(M), the assertion is evident. Furthermore, if A = 72^2 7^ 0, then 

% )fc) $ 2 (A ,0 m,h)} = (72 — 7l )/ + /C/i for all (/,/») G (C e 3 +“ (M)) 2 , 

whereby 

K,y^b k cos(R k x) =-AV . . f* . . b k cos(R k x). 

^ sinhUti-fflo 
fceN fceN 1 * ' 
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It is not difficult to see that /C(Cg+^.(R)) C C^ er (M), so that Imc^^^lA, 0) = Cgj^(M). 
Hence, Im9y ^$(A, 0) is not a closed subspace of (Cgp^[.(M)) 2 . This proves (z). 

To prove (zz), we choose A 0 { 0 , 72 ^ 2 } and set 

D(k, A) := ml 1 (A)ml 2 (A) - mf (A)mf (A), k <E N. 

From (2.13)-(2.16) it is clear that there exists ko € N such that D(k, A) 7 ^ 0 for all k > ko- 
Defining the symbols m 1J (A) by mV (A) = mV (A) for k > ko and 1 < i, j < 2 and 

mj^A) = m| 2 (A) = 1, m^, 2 (A) = m^ 1 (A) = 0 for 0 < k < ko — 1, 

we see that 0 ) is a compact perturbation of the operator 

T ■= ~^ 2 ) € £((C e 3 +“ (M)) 2 , (Cg 2 +“ (M)) 2 ), 


where Ajj € £(C 3 +“ r (M), C 2 +" r (M)) is the Fourier multiplier defined by rrd J (A), 1 < z, j < 2. 
Because D(k,A) := fhj}(A)fh 22 (A) — m^. 2 (A)m| 1 (A) / 0 for all k € N, we can define the 
formal inverse of T by 

^ / -Bn Bi2\ 

\B 2 i B 22 J 

Here, B^ is the Fourier multiplier corresponding to the symbol b 13 , 1 < i,j < 2, whereby 


b 


11 

k 


g 2 (A) b i2 . g 2 (A) . mg 1 (A) ^2. ffl(A) 

D(k,A)' k ' D(k, A)’ ' D(k,A)’ k ' D(k, A) 


for k £ N. 


Using now [21, Theorem 2.1], we see that a Fourier multiplier 


y: a k COS (R k x) -A- AfcOffc cos (B fc x) 
fceN fceN 


belongs to C (C 2 ^ r (K), Cf^r W) if 

sup | A:Afe | < 00 and 
km 


sup k 2 \X k+1 - Afe| < 00 . 

fcgN 


Because of this, it is a matter of direct computation to see that Bij € £(C 2 p^,(M), C 3 p^.(M)) 
for all 1 < i,j < 2. Hence, T is an isomorphism, and therefore d^fj j)<k(A, 0) is a Fredholm 
operator of index zero. □ 


Because of Lemma 2.5 (i) it is clear that we cannot apply the Crandall-Rabinowitz bi¬ 
furcation theorem at (A, 0) with A € { 0 , 72 ^ 2 }- As a consequence of this, the laminar flows 
from which we show that non laminar waves bifurcate will not possess stagnation points 
at the wave surface or on the interface separating the two layers of constant vorticity, cf. 
(2.6)-(2.7), but only inside the layers. This is different than in the case of internal waves 
propagating between two layers of constant but different density, where in the presence of 
capillarity stagnation points may be located also on the internal wave, cf. [27]. 

It is now evident that potential bifurcation values for A 0 { 0 , 72 ^ 2 } are to be looked for 
among the solutions of 

D(k, A) = 0 (2.17) 

for some integer k > 1. Since in Theorem 2.3 the Frechet derivative d(f ^^(A, 0) needs to 
be a Fredholm operator of index zero with a one-dimensional kernel, we need to find A such 
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that (2.17) has exactly one root 1 < k G N. Plugging the expressions (2.13)-(2.16) in (2.17), 
we rediscover the dispersion relation 


A 3 — -i- 


72 


Rk 

+ tanh (R k d) 
tanh (R k d) 


smh(R k d 2 )cosh(R k di)\ 
Rkd2 + cosh (R k d) ) 


+ 7i 


sinh(R k di) cosh(R k d 2 ) 


l2 d 2 ~ 9 

Rk 


+ 72(71 -72) 


cosh (R k d) 

sinh(R. k di ) sinh(i?*.d2)" 


A 2 


R'i sinh (R k d) 


+ 9- 




(7i - 72 ) 


sinh(R, k d\ ) sinh(R k d 2 ) 


+ r )2 d 2 Rk 


A 


= 0, 


(2.18) 


sinh (R k d) 

found also in [25, Equation (5.11)] (with a = 0). This relation has been analyzed in the 
setting of flows without stagnation points in [5] for 71 7^ 0 = 72 and in [9] for 71 = 0 7^ 72. 
Herein, we assume only that 71 7^ 72 and restrict the analysis to the complementary case 
when stagnation points are included. 

I11 studying the dispersion relation (2.18) we will make use of the following remark, which 
allows us to restrict our attention to a few of relevant cases, the remaining ones being 
analogous. 


Remark 2.6. Note that (2.18) possesses the following symmetry property: A; is a solution 
of (2.18) for some A 0 {0,72^2} and (71,72) € R 2 if and only if A; is a solution of (2.18) 
for —A 0 {0, —72^2} and (-71,-72) € R 2 . Because additionally the inequalities (2.6) and 
(2.7) are invariant under the transformation (A, (71,72)) H > (—A, (—71,—72)), we are left 
only with the two cases: 

(*) 72 > 0 and 71 / 72; 

(ii) 72 = 0 and 71 < 0. 


3. Analysis of the dispersion relation: the case 72 > 0 and 71 7 ^ 72 


Because the dispersion relation is highly nonlinear in k, the study of the roots of (2.18) 
when keeping A fixed seems to be very difficult. Therefore, we consider the inverse problem 
of determining the zeros Ai, A2, A3 of this cubic equation when keeping k fixed, and then 
to study the properties of the mappings k i->- Aj(A;), 1 < i < 3. We will do this for small 
wavelength L, because then we can use asymptotic expansions and Cardano’s formula in 
order to determine the roots A* of (2.18). This small wavelength regime corresponds to the 
setting t —>• 00, where 


t := R k 


2-irk 

— 


€ R. 


Plugging t. for R kl the equation (2.18) can be written in the more concise form 


A 3 + A(t) A 2 + B(t) A + C(t) = 0. 


(3.1) 


We will show in the sequel that equation (3.1) has three real roots when t is sufficiently 
large. To this end, we note first that the coefficient functions A = A(t), B = B(t), and 
C = C(t) and their first derivatives have the following asymptotic expansions for t —>• 00: 


A = -72CA2 - 7172 +o(^T 


2t 


t 2 


D ll d 2-9 , 72(71-72) , /l 

B = -;- 1 -—75-h O -r 


t 


2 1 2 




_ 57 2 d 2 5(71 - 72) . /I 

C — ;-1-—^-r o 


2t 2 


t 3 


A' = 


71+72 1 


B' = — 
C' = - 


2 t 2 
l\d 2 - g 




t 2 

572^2 

t 2 


72(71-72) , /I 

- + O TT 


t 3 

5(71 - 72) 

t 3 


t 4 


+ o 


(3.2) 
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Letting z := A + .A/3, we find that 2 solves the depressed cubic equation 

„3 


with 


P 

3 


B 

3" 


A* 

9 


z + pz + q = 0 , 

(72^2) 2 , 72^2(272 - 71) - 3g 47172 


(3.3) 


9 


+ 


9t 


+ 


m 2 


+ c| ¥ 


and 


q _ A 3 AB C 

2 _ 27 _ IT + "2 

( 72 ^ 2) 3 


27 


+ 72^2 


6g +72^2(272 - 71) 
187 


' , 7i ~ 7172 + 72 , g(272 -7i) l J_ 

72 2 36 + 3 J t 2 


| 972(7^ - 7I) - (71 + 72) 3 1 | 3 / 1 


216 7 3 

Observe that the discriminant for (3.3) is 


f 4 


D : = 


P 


+ 


9g (72 cfc)' 4 
2437 


+ o( 


<0 for 7 —> 00 , 


property which implies, cf. [33], that (3.3), and hence also (3.1), has three real roots. They 
are given by the relation 2 = r cos(/3), whereby 


r = 


-4 p 272^2 72^2(272 - 71) - 3g 


and j3 is one of the solution of 


cos(3/3) = --\ -o = 1 - 


27 


372^2^ 
3 3 <? 


+ 0 


(?) 


(3.4) 


+ 0 


2 ( 72 ( 72) 2 7 ^ 7 2 


(?)• 


Thus, choosing /3 := 3 1 arccos ( — (g/2)-^/—27/p 3 ) we see that /3(7) \t _ > . 00 0 and the roots 
of (3.1) are 

^4 

Ai = r cos(/3) — 


Aa = 


r cos 


+ 


3’ 



27T\ 

A / 


— 

-= —r cos 

ft + - 

3 7 

3 V 

P 3 

2 vr\ 

A ( 

„ 7T 


3 7 3 ' 


.1 

3"’ 

A 


(3.5) 


Together with (3.2) and (3.4), it follows at once that for t —> 00 we have 

A-i(t) —>■ 72(72, ^ 2 ( 7 ) —> 0, Aa(t) —>• 0. (3-6) 

Let us also observe that since C{t) > 0 for t —>• 00 , it must hold that A. 2 (t)A${t) < 0 for 
7 -A 00 . Moreover, it is clear from (3.5) that A 2 ( 7 ) > A 3 (f), hence 

A 2 ( 7 ) > 0 > A 3 (7) for t —> 00 . 

3.1. Existence of water flows bifurcating from Ai. In order to consider the bifurcation 
problem for ( 2 . 11 ), we need to study first the properties of the mapping [t e-)■ Ai(t)]. 

Lemma 3.1. There is a constant to > 0 such that the function 

[[fo, 00 ) 7 t eA Ai(t) € (0, 00 )] 


is strictly monotone. 
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Proof. Let f>(t, A) := A 3 + Aft) A 2 + B(t )A + C(t) for A € R and t > 0. Since for t —t oo we 
have 

0 a(L Ai(f)) = (Ai(t) - A 2 (t))(Ai(t) - A 3 (t)) > 0, (3.7) 

we conclude that Ai is differentiable with respect to t. On the other hand 


t 2 $ t (t,Ai(t)) = t 2 (A 2 1 (t)A'(t) + A 1 (t)B'(t) + C\t)) -+ t _>oo 


(72^2) 2 (7i - 72 ) 

2 


Since A^(t) = — <j>t(t, Ai(t))/<^>A(t, Ai(t)), we conclude that A( has the same sign as 72 — 71 - 
The constant to is defined as to := inf{£ > 0 : |A^| > 0 on (t, 00 )}. □ 


From Lemma 3.1 it follows at once that 

• if 71 < 72 , then Ai(t) satisfies (2.6) for t > to; 

• if 71 > 72 , then Ai(t) satisfies (2.7) for t > to- 

We look now for bifurcation solutions when choosing Ai as the bifurcation point. There¬ 
fore, we choose to > 0 in Lemma 3.1 large enough to guarantee additionally that 

inf A 2 > sup (A 2 + A 3 ), 

[to ,o°) [to, 00) (3.g) 

17(0, Ai(t)) / 0 for all t > to- 
Let 

L 0 := 27t/£ 0 , (3.9) 

fix L < Lq, and set Ai := Ai(2ir/L). Then, since cj)(2ir/L, A\) = 0, we get D(l, Ai) = 0. 
Due to the choice of to, the equation D(-, Ai) = 0 has no solutions k £ N other than k = 1. 
Consequently, since Ai ^ {0, 72 ^ 2 }? the derivative d^f t h)^i(A\, 0) is a Fredholm operator 
with a one-dimensional kernel 

Ker 0(j- / l )<I > i(Ai, 0) = span {(m 22 (Ai), — m 21 (Ai)) cos(27rx/L)|. (3.10) 

In order to apply Theorem 2.3 to this particular setting, it remains to study whether the 
transversality condition is satisfied. To this end, we obtain the following characterization of 

i(Ai,°). 

Lemma 3.2. Let Lq be given by (3.9), L < Lq, and set A\ : = Ai(2n/L). Then, we have 
Im %ft) $ ( A i>°) = = (£a cos (M,^ % cos(M) : 7i = 

(3.11) 

Proof. To prove the claim, let (/, h) = (X)fceN ak cos i^kx ), XlfceN cos (Rkx)) be such that 
% ift) $(Ai,0)(/,/i) = (£, 77 ). Then, obviously 

( ml 1 (Ai)ai + m\ 2 (A 1 )bi = 'y 1 , 

\ ml 1 {A 1 )a 1 + mf(Ai)bi = tji. 

Because D(l, Ai) = 0, we find from (2.13)-(2.16) that 

mfjAi) _ mP(Ai) _ 
m 21 (Ai) mf 2 (Ai) 

Hence, (£, 7 ) is an element of the set defined by the right-hand side of (3.11). Because the lat¬ 
ter set is a closed subspace of (C'gj^-(M)) of codimension one that contains Im9(j j / l )$(Ai, 0 ), 
the conclusion follows from Lemma 2.5. □ 


We are now at the point of showing the transversality condition (c) from Theorem 2.3. 
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Lemma 3.3. We have that 

5A(/,h)^(A 1 ,0)[(mf 2 (Ai),—mf(Ai)) cos(2vrx/L)] 0 Imd(/,fc)$(Ai,0). 

Proof. Since for a, b € M, it holds that 

5a(/,/i)^(Ai, 0) [(a, 6) cos(2vrx/L)] 

= (am}J A (Ai) + bm\ 2 A (Ai), am^ A (Ai) + 6 m 22 A (Ai)) cos(2irx/L ), 
we are left to show that 

mf (Ai)m^ A (Ai) -mf(Ai)mJ 2 A (Ai) (mf (Ai)m^ A (Ai) - (Ai)m^ 2 A (Ai)), 


or equivalently that 

m 22 (Ai)m^ A (Ai) - m 21 (Ai)mi 2 A (Ai) / mJ 2 (Ai)m 2 ^ A (Ai) - mJ 1 (Ai)m 22 A (Ai). 

Hence, we need to show that D\(l, Ai) 0. Recalling the definition of the mapping from 
the proof of Lemma 3.1, we have that 11(1, A) = /L, A), and therefore 


£>a(1, Ai) = $ A (27r/L,Ai(27r/L)) > 0, 


which is the desired property. □ 

Theorem 3.4 (Bifurcation from Ai). Let 72 > 0, 71 7 ^ 72 and let a € (0,1) be given. 
Furthermore, let Lq be the constant defined by (3.9) and L < Lq. Then, there exists a real- 
analytic curve (A, (/, h)) : (—e, e) —>• (0, 00 ) x O consisting only of solutions of the problem 
(2.11). This curve contains exactly one trivial solution of (2.11), and for s —>• 0 we have 
that 


A(s) = Ai + 0(s), (/, h)(s) = s(m 22 (Ai), -m 21 (Ai)) cos(27r x/L) + 0(s 2 ), 

whereby Ai := Ai(27 t/L). The flow determined by (A (s),(f,h)(s)),s € (—e,e), contains a 
critical layer consisting of closed streamlines very close to the internal wave 

(■ i ) in the layer adjacent to the wave surface if 71 <72, or 
(ii) in the bottom layer i /71 > 72 . 

Moreover, the amplitude of the internal wave is much larger than that of the surface wave, 
cf. Figure 1. 


Proof. It remains only to show that the amplitude of the internal wave is much larger than 
that of the surface wave. To this end, we note that due to D(l, Ai) = 0, we have 


since 


m{ 2 (Ai) _ m\ 2 { A t ) 

mf 1 (Ai) mp(Ai) 


sign ( 71 - 72 ) 00 , 


m 


L-Somi'lAi) 


(A,) _ g + ^W-A 
- Al (()( 72 *-A 


= lim 
t —>00 


= sign (72 - 71 ) 00 . 


□ 
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FIGURE 1. This figure illustrates the streamlines in the moving frame for the 
solutions found in Theorem 3.4 for 71 > 72 > 0 (left) and 71 < 72,72 >0 (right), 
cf. Lemmas A.1-A.2. The thick streamlines represent the wave surface, the internal 
wave, and the flat bad, respectively. The blue streamlines are separatrices which 
bound the critical layer and the dashed line consists of points where the y —derivative 
of the stream function vanishes. This line contains in both cases exactly three 
stagnation points: two located at x = 0 and x = L, and a third one inside the 
critical layer at x = L/ 2 . 


3.2. Existence of water flows bifurcating from A 2 . For t —>• 00 we have that 

• A 2 (t) satisfies ( 2 . 6 ); 

• if 71 di + 72^2 £ 0, then A 2 (t) satisfies also (2.7). 

Letting <p = (f)(t, A) be the function defined in the proof of Lemma 3.1, we note that for 
large t we have 

0 a(L A2(t)) = (A 2 (t) — Ai(t))(A2(f) — A3(f)) < 0. 

Hence, A 2 is differentiable with respect to t. Moreover, it follows from (3.2) that 

f 2 0 t (t, A 2 (t)) -H-kxj -£72^2, 

and therefore A ' 2 (t) = — A 2 (t))/</>A(t, ^(t)) < 0 when t is large. Defining 

to := inf{t >0 : A' 2 < 0 on (t, 00 )}, 

we see that [[to, 00 ) A 2 (f) € (0, 00 )] is decreasing. In view of this property, we can 

choose to > 0 large enough to guarantee that 

sup A 2 < inf Ai, 

[to,00) [A, 00 ) ^2 12^ 

D(0, A 2 (t)) / 0 for all t > to- 


Then, we set 


Lq 2n/to, 


(3.13) 


we fix L < Lo, and define A 2 := A2(27r/L). Recalling that (j>(2n/L, A 2 ) = 0, we obtain that 
D(1,A2) = 0. In fact, the equation D(-, A 2 ) = 0 has k = 1 as the only integer solution, 
cf. (3.12). Because of A 2 € ( 0 , 72 ^ 2 ), i(A 2,0) is a Fredholm operator with a one¬ 

dimensional kernel 


Ker 9 (j fe )$i(A 2 ,0) = span {(m 22 (A2), — m 21 (A 2 )) cos( 27 rx/L)}. 














14 


C. I. MARTIN AND B.-V. MATIOC 


Using the same arguments as in the proof of Lemma 3.2, we see that 


Im %A)<L( A 2,0) = | (£,77) = 6c cos {Rkx), ^2 rj k cos(i4x)) 

keN ke N 

Moreover, the transversality condition 


: 71 = 


m} 1 (A 2 ) I 


0A(M)^(A 2 ,O)[(mf (A 2 ),-m? 1 (A 2 )) cos( 2 vrx/L)] 0 Im9 (/>fe) $(A 2 ,0) 

reduces to showing that L>a( 1 ,A 2 ) = c/)\(2it/L, A 2 ( 27 t/L)) 7 ^ 0 , relation which holds true. 
We conclude with the following result. 




FIGURE 2. This figure illustrates the streamlines in the moving frame for the 
solutions found in Theorem 3.5 for 71^1 + 72^2 > 0 (left) and 71 di + 72^2 < 0 
(right), cf. Lemmas A.3-A.4. 


Theorem 3.5 (Bifurcation from A 2 ). Let 7 2 > 0, 71 7 ^ 7 2 and let a € (0, 1) be given. 
Furthermore, let Lq be the constant defined by (3.13) and L < Lq. Then, there exists a real- 
analytic curve (A, ( f,h )) : (— e,e) —» (0, 00 ) X O consisting only of solutions of the problem 
( 2 . 11 ). This curve contains exactly one trivial solution of ( 2 . 11 ), and for s —>■ 0 we have 
that 

A(s) = A 2 + 0(s), (/, h)(s) = s(mf(A 2 ), -m 21 (A 2 )) cos(2nx/L ) + 0(s 2 ), 

whereby A 2 := A 2 (27r/L). The flow determined by (A(s), (/, h)(s)), s € (— £,e), contains a 
critical layer consisting of closed streamlines 

(i) in the layer adjacent to the wave surface 2 / 71^1 + 7 2 d 2 > 0; 

(ii) in each of the layers 1 / 71^1 + 7 2 d 2 < 0 . 

The vortex in the top layer is located right beneath the wave surface. Moreover, the amplitude 
of the internal wave between the two layers is much smaller than that of the surface wave, 
cf. Figure 2. 


Proof. It remains only to show that the amplitude of the surface wave is much larger than 
that of the internal wave. This follows from (2.15)-(2.16) , as we have 


A 2 ) 

mf{ A 2 ) 


— lL->0 —OO. 


□ 
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3.3. Existence of water flows bifurcating from A 3 . Since 0 > A 3 (f) -A t~>oo 0 we see 
that A 3 (t) satisfies (2.7) provided that 71^1 + 72^2 < 0. Because for large t 

<M*,A 3 (f)) = (A 3 (t) - Ai(t))(A 3 (t) - A 2 (t)) > 0, 

the function A 3 is differentiable with respect to t. Since t 2 (j)t(t, A^(t)) —tt-too ~ 572 ^ 2 , we 
conclude that A ' 3 (t) > 0 when t is large. Defining to := inf{t >0 : A 3 > 0 on (t, 00)}, we 
see that [[to, 00) 9 i O A ^(t) € (—00, 0)] is increasing. In view of this property, we can 
choose to > 0 large enough to guarantee that 

.0(0, A^ft)) 0 for all t > to- (3-14) 

Let 

L 0 := 27r/t 0 , (3.15) 

choose L < Lq , and define A 3 := /L). Since c))(2tt/L, A 3 ) = 0, we get that k € N solves 

D(k, A 3 ) = 0 if and only if k = 1. Moreover, since A 3 < 0, ^1 (A 3 ,0) is a Fredholm 

operator with a one-dimensional kernel 

Ker<9 (/ft) $i(A 3 ,0) = span {(mf (A 3 ), -m 21 (A 3 )) cos(2t tx/L)}. 

As in Lemma 3.2, we find that 

Im %,/i)4>(A 3 ,0) = {(£, 7 ) = (^2^ k cos(R k x),^2r] k cos(R k x)j : 71 = 

keN feeN mi ^ 3 ' 

the transversality condition 

^A(/,h)4 > (A 3 , 0 ) [(m 22 (A 3 ), —m 21 (A 3 )) cos(2vrx/L)] ^ Im %,fc)$(A 3 ,0) 

being equivalent to D\(l, A 3 ) = /L, A^(2tt/L)) 0. This shows that all the assump¬ 

tions of Theorem 2.3 are satisfied. Consequently, we have the following result. 

Theorem 3.6 (Bifurcation from A 3 ). Let 72 > 0, a € (0,1), and assume d\ + 72^2 < 0. 
Furthermore, let Lq be the constant defined by (3.15) and L < Lq. Then, there exists a real- 
analytic curve (A, (f,h)) : (— e,e) (0, 00) x O consisting only of solutions of the problem 
(2.11). This curve contains exactly one trivial solution of (2.11), and for s — >• 0 we have 
that 

A (s) = A 3 + O(s), (/, h)(s) = s(m 22 (A 3 ), -m 21 (A 3 )) cos( 2 vr x/L) + 0(s 2 ), 

whereby A 3 := A^{2-k/L). The flow determined by (A (s),(f,h)(s)),s € (— e,e), contains a 
critical layer consisting of closed streamlines in the layer adjacent to the bed. Moreover, 
the amplitude of the internal wave between the two layers is much smaller than that of the 
surface wave, cf. Figure 3. 


Proof. The claim concerning the amplitude of the surface and internal waves follows from 
(2.15)-(2.16), as we have 


mf{ A 2 ) 
mf (A 2 ) 


>0 00 • 


□ 
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FIGURE 3. This figure illustrates the streamlines in the moving frame for the 
solutions found in Theorem 3.6 (left) and Theorems 4.2 and 4.3 (right), cf. Lemmas 


A.5-A.6. 


4. Analysis of the dispersion relation: the case 72 = 0 and 71 < 0 

Because of 72 = 0, the inequality (2.6) reduces to A = 0, situation when 9(/ i / l ) < h(A, 0) is 
not even a Fredholm operator, cf. Lemma 2.5. For this reason we consider the bifurcation 
problem for (2.11) just for values of A which satisfy (2.7). Hence, the flows that we construct 
will have stagnation points only in the bottom layer. 

With the notation from Section 3, we determine for the depressed cubic equation (3.3) 
that 

D = (|) + (|) = - f ^” 3 + °i t ~ 4 ) < 0 for t —» 00 , 

hence (3.1) has again three positive roots. They are given by the relation 2 = rcos(/3), 
whereby 


r =y?= 2 #" /2+ i^ 3/2+ °T 5/2 ) <“> 

and /3 is one of the solutions of 

- -iff - -^r 12 +ff - 3/2 + 

Setting /3 := 3 _1 arccos ( — (g/2)y / —27/p 3 ), we see that (3(t) 00 7 t/6 and the roots of 

(3.1) are 

A 

Ai = rcos(/3) - —, 

A 2 = rcos (p- y )-j = —r cos (/? + |) - j, (4.2) 

A 3 = rcos (<3 + f )-f = - rcOS (/3 - f ) - f. 
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It follows now easily from (3.2), (4.1), and (4.2) that Aj 0 for i € {1,2,3} and that 

A 3 < A 2 < 0 < Ai for t —>• 00 . Thus, we can find to > 0 such that 

7 idi < A 3 < A 2 < 0 < Ai on [t 0 , 00 ), 

14(0, Aj(t)) 7 ^ 0 for all t > to, i = 2, 3. 

I 11 view of 71 < 0, the relation (2.7) is equivalent to A € ( 71 ^ 1 ,0) and therefore just 
flows bifurcating from negative A may contain stagnation points. For this reason, we only 
investigate in the following the functions A 2 and A 3 . 

Lemma 4.1. There exists to > 0 such that (4.3) holds and Aj : [to, 00 ) —>• (— oo,0),i = 2,3, 
are both increasing functions. 


Proof. Note first that (f\(t, A 2 (t)) < 0 and <f\(t, A^(t)) > 0 for t > to- Hence, Aj,i = 2,3, 
are differentiable on [to, 00 ). Moreover, it is easy to see from (4.1) and (3.2) that 

lim t 5 / 2 <Mf,A 3 (f)) = Hm t 5 / 2 A 3 (t)H'(t) = -g^ 2 . 

£—>•00 t —>00 

Hence, we may chose to large to ensure the assertion for the mapping A 3 . 

This argument does no longer work for A 2 as cos(/3 + 7 r/ 3 ) —>•*_>.oo 0. Hence, we have to 
determine an expansion for cos (ft + 7 r/ 3 ). Let 


and observe that 



7i 

6 



| cos(/3) — -^0 11 cos 2 (/3) + zq cos(fl) + Zq\ =| 4 cos 3 (/ 3 ) — 3 cos (/3) — 4 zq + 3 ^o I 

cos(3/3) + 71 {/“t _1 ^ 2 + 0(t _1 ) = 0(t _1 ). 


Therewith cos(/3) = zq + 0{t 1 ) for t — >• 00. It is now easy to see that for t —>• 00 we have 

sin(/3)=i + -^t- 1 /2 + 0 (t-i ) , 

2 2 y/g 

from which it follows easily that 

A 2 (t) := “^t -1 + 0(t~ 3 / 2 ) and A 3 (t) := ~^t~ 1/2 - ^t~ x + 0(t~ 3 / 2 ). (4.4) 

fa fa 

The expansion (4.4) combined with (3.2) shows that 

lim t 3 <ft(t, A 2 (t)) = lim t 3 (A 2 (t)H'(t) + C'(t)) = —371/6 > 0 , 


relation which proves the claim. 


□ 


Theorem 4.2 (Bifurcation from A 3 ). Let 72 = 0, 71 < 0, and a € (0,1). Furthermore, 
let Lq \= 2ir/to and L < Lq. Then, there exists a real-analytic curve (A, ( f,h )) : (— e,e) — » 
(0,oo) x O consisting only of solutions of the problem (2.11). This curve contains exactly 
one trivial solution of (2.11), and for s —>• 0 we have that 

A (s) = A 3 + O(s), (/, h)(s) = s(m 22 (A 3 ), -m 21 (A 3 )) cos(27 tx/L) + 0 (s 2 ), 

whereby A 3 := A 3 (27r/L). The flow determined by (A(s),(f,h)(s)),s € (—e,e), contains a 
critical layer consisting of closed streamlines in the layer adjacent to the bed just below the 
internal wave. Moreover, the amplitude of the internal wave between the two layers is much 
smaller than that of the surface wave, cf. Figure 3. 
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Proof. Because t : = 27r/L > to, we know from (4.3) and Lemma 2.5 that $(/,/i)^i(A 3 ,0) is 
a Fredholm operator. To determine its kernel we need to solve D(k, A 3 ) = cf>(kt, A3) = 0. 
As A3 = A3 (t), we see that D(1,A3) = 0, while (4.3) ensures that D( 0, A3) = 0. Recalling 
Lemma 4.1, we see that D(k, A3) f 0 for all k > 2. Indeed, for k > 2, A^(kt) > A3, and if 
A3 = A2 (kt), then A3 = A2 (kt) > A2 (t) > A3, a contradiction. Hence, d^f ^^i(A^, 0) is a 
Fredholm operator with a one-dimensional kernel 

Ker < 9 (/,fe)$1(A 3 , 0 ) = span {(mf 2 (A 3 ), -mf 1 (A 3 )) cos(2vrx/L)}. 

Similarly as before we have 

lmd (f,h)$(A 3 , 0 ) = {(£, 77 ) = (^ 2 ^ k cos(R k x),^ 2 r] k cos(R k x)j : 71 = 

fceN feeN mi ^ 3 ' 

and one can verify that the transversality condition 

0 A(M)^(A 3 ,O)[(mf (A 3 ),-m? 1 (A 3 )) cos(2ttx/L)] 0 Im 9 (/>fe) $(A 3 , 0 ) 


is also satisfied. We are thus in the position of applying Theorem 2.3. To finish the proof, 
we infer from (2.15)-(2.16) and (4.4), that 


wp(A 3 ) 

mf (A 3 ) 


~tL-s-0 00 • 


□ 


When considering bifurcation from A 2 the situation is more complicated because the 
derivative 3 (j/ i ) < hi(A 2 ,0) may possess a two-dimensional kernel if A%(2t xk/L) = A2(2ir/L) 
for some L > Lq and some integer k > 2. When this happens, the integer k is unique, cf. 
Lemma 4.1, so that we can conclude the existence of a curve of bifurcating solutions from 
Theorem 4.2. When As(2nk/L) f A2(2 tx/L) for all k > 2, we can apply again Theorem 2.3. 

Theorem 4.3 (Bifurcation from A 2 ). Let 72 = 0, 71 < 0, and let a € (0,1). Furthermore, 
let Lq := 2n/to and L < Lq. 

(i) Assume that A^,{2 t: k/L) 7 ^ A2(2tt/L) for all integers k > 2. Then, there exists a 
real-analytic curve (A, (/, h)) : (—e, e) —>• (0, 00 ) X O consisting only of solutions of 
the problem (2.11). This curve contains exactly one trivial solution of (2.11), and 
for s —} 0 we have that 

A(s) = A 2 + O(s), (/, h)(s) = s(mf(A 2 ), -mi 1 (A 2 )) cos(27 tx/L) + 0(s 2 ), 

whereby A 2 := A 2 ( 2 n/L). 

( ii ) Assume that Aj,(2i{k/L) = A 2 for some integer k > 2. Then the assertion of Theorem 
4-2 holds true, but with L replaced by L/k. 

The flow determined by (A (s),(f,h){s)),s € (—e,e), contains a critical layer consisting of 
closed streamlines in the layer adjacent to the bed just beneath the internal wave. Moreover, 
the amplitude of the internal wave between the two layers is much smaller than that of the 
surface wave, cf. Figure 3. 

Proof. Setting t := 2 -k/L > to, we know from (4.3) and Lemma 2.5 that 3 (j j /j)<I > i(A 2 , 0) is 
a Fredholm operator. To determine its kernel we need to solve D(k, Af) = &{kt, A 2 ) = 0. 
A solution of this equation is k = 1 as A 2 = A 2 (t). Equation (4.3) ensures additionally that 
0(0, A 2 ) = 0. Because A 3 is increasing to zero, there may exist a (unique) integer k > 2 
such that As(2ixk/L) = A2(2tt/L), hence D(k, A 2 ) = 0. In this case we are in the situation 
(ii) and the proof is obvious. If Az(2irk / L) f A2(2 it/L) for all integers k > 2, then we are 
in the case (i) and the proof is similar to that of Theorem 4.2. □ 
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Remark 4.4. Since the properties of the functions Aj(f), for i = 1,2,3 were essential in 
finding the branches of solutions to the water wave problem, we summarize them in the 
Table 1 below. 


- 

7i > 0 

71 = 0 

7i < 0 

72 > 0 

A 3 < 0 < A 2 < Ai 

Ai —> r y 2 d 2 , Aj —> 0, i £ {2,3} 

0 

II 

CM 

O 

A 3 < 0 < A 2 < Ai < 7 idi 

Aj -a- 0, *£{1,2,3} 

- 

71 di < A 3 < A 2 < 0 < Ai 

A* —> 0 , *£{ 1 , 2 ,3} 

72 < 0 

A 3 < A 2 < 0 < Ai 

A 3 — > r j 2 d 2 , A* —» 0 , i € { 1 , 2 } 


TABLE 1. Properties of the roots Aj,i £ {1,2,3}, of the dispersion relation 
(2.18) in dependence of the vorticity constants 7*,? £ {1,2} for 71 7^ 72 and 
for large R.f- = 2nk/L. Our analysis is dedicated to the cases: ( i ) 72 > 0 and 
71 7^ 72; and (ii) 72 = 0 and 71 < 0. The analysis in the other two cases: 
(***) 72 < 0 and 71 7^ 72; and ( iv ) 72 = 0 and 71 > 0 is similar to that for 
(i) and (ii), respectively (see Remark 2.6). 


Appendix A. 


We present herein the proof of Lemma 2.4 and additionally we rigorously prove that the 
streamline pattern for the solutions that we found is as shown in Figures 1-3, respectively. 
To this end, we determine first explicit expressions for the elliptic and boundary operators 
introduced right before Observation 2.1. Given (/, h) £ O, it is easy to see that 

Mf)~d 2 d + V f'd | 1 (d + y)f'\ 2 + d\ (d l + f)f"-2f 2 

A{f)-a xx i di + f S0 xy + ( di + f ) 2 cfo W + y) ( di + /)2 


■'yi 


(A.l) 


A(f, h) = d xx - 2 


d 2 ti + (h-f)'y \d 2 h! + (h - f)'y\ 2 + d% 


h- f + d 2 
d 2 h" T (h f)"y 


dxy T 


(h-f + d 2 f 


d, 


yy 


- 2 


d 2 h'(h - /)' + (h' - ffy^ 


d y , 


h-f + d 2 ~ (h-f + d 2 y 

respectively, given (w\,w 2 ) € C'^f° r ({},\) x C^^ r (Ll 2 ) and A £ R, we have that 


(A.2) 


Bi(A,(f,h),w 2 )= \d x w 2 \ 2 - 


#2 (A, (f,h))[wi,w 2 ] = 

Proof of Lemma 2-4- Since 


h- f + d 2 
+ 2 g(d + h) — Q( A), 
d 2 


2d 2 h! d 2 2 (l + h' 2 ) 2 

d x w 2 d y w 2 + —- , 2 \dy w 2\ 


(h- f + d 2 ) 2 ' 


y =0 


h - f + d 2 


d y w 2 


di 


y=-d ,2 d\ + / 


d y W 1 


y=-d 2 


(A.3) 

(A.4) 


,0) [(/,/*)] 


d f MA , 0) [/] d h ^(A ,0)(h] 
d f * 2 (A , 0) [/] 9 h 4» 2 (A,0)[h] 
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we only need to determine the entries in the matrix 0). 

The derivative <9/<hi(A,0) Using the definition of d>i, we see that 


/ 


<9/d>i(A, 0)[/] =2 — \d y ^ 2 \- + d y ^dy(d f w 2 (A,0)[f}) 
1 «2 


y =o 


(A.5) 


whereby | y=0 = A and z := dfW 2 {k, 0)[/] is, in view of (2.10), the solution of the 
Dirichlet problem 

Az = — <9/71(0,0)[/]-02 in U 2 , 
z = 0 on <9D 2 - 

A routine calculation shows now that 


(A.6) 


3/71(0,0) [/]^ = 


272 / 

do. 


72 2 


A 


+ KV + ir// / • 


do 


Expanding / and z(y),y G [—c? 2 ■> 0], by their Fourier series, we have 

f = ^2 a k cos (R k x) and z(y) = ^ a k A k (y) cos (R k x). 
he N fceN 

The coefficients A k solve, in view of (A.6), the following boundary value problem 

/ K-RlAk=-% + + T,y) in (-d 2 ,0). 


and therefore 


A k (—d 2 ) = Afc(0) = 0, 

Mv) = (A - 7 ! d 2 )4 nh(fla ' ) 


72 2, A 

V + V 


sinh(i?fcC?2) 

Using the relation (A.5) we obtain now that 

a/*i(A,o)tfl = E m]}a k cos {R k x), 


ke N 

whereby (m). 1 )fceN is defined by (2.13). 

The derivative 0) We have that 

<9fo$i(A,0)[/i] = 2 — -^-1^21 2 + dy^2dy{dhW 2 {A,0)[h\) 
l a 2 

with z := dh.w 2 (A, 0)[/i]) solving the Dirichlet problem 

Az =—0^71(0, 0)[/i]V’2 in D 2 , 

z = 0 on (AI 2 , 

cf. (2.10). Recalling (A.2), we compute that 


y =0 


+ 2 . 9 / 1 , 


3^71(0, 0)[/i]^2 {^ y2 + + A ) h 


72 2 , 72^2 + A 


d 2 - 


d 2 


Using Fourier expansions as before, that is 

h = b k cos (R k x) and z(y) = ^ h k B k (y) cos (R k x) for y G [~d 2 , 0], 
fceN feeN 

we obtain that the coefficients B k satisfy 

f B'l - R\B k = %- Rl(jy + + A ) in (-*,0), 

1 B k (—d 2 ) = B k { 0) = 0. 
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The solution of this boundary value problem is 

/ sinh (R k y) \ j 2 2 . 72^2 + A 

Bt{y) = - A ( tanh (i4d 2 ) + C ° Sh {Rty) ) + di y + ^^ y + A ’ 
and the desired representation for the derivative 0) follows at once. 

The derivative df& 2 (A, 0) From the definition of $2 we obtain that 

d/$ 2 (A,0)[/] = \^-d y ip% + J t-dy^ + dy(dfW 2 (A,0)[f]-d f uj 1 (A,0)[f])] 

L«2 «1 j y=—d2 

whereby, in the equality above, 2 := dfW\(A,0)[f] solves the Dirichlet problem 

A z = — dfA(0) [f]ipi in fii, 

on dQ \. 


z = 0 

In view of (A.l), we compute that 


dfAmm=- 


271 / 

d\ 


fy 2+ 

d\ 


d + d 2 A- 7 2 d 2 \ dA dd 2 (7i— 72 ) 


di 


-71 + 


di 


) y+ i7 + 


d 1 


Expanding / and z(y),y € [—d, — d 2 \, by their Fourier series 

/ = ^2°k cos (R k x) and z(y) = ^ a k C k (y) cos (R k x) 
fceN feeN 

we find that the coefficient C k is the solution of 

i/l r>2 


di k Gfc — 


— ^71 _ p2 
di n k 


TL y 2 _|_ 
di» ^ 


( d+d2 I A-72d2 

l di 71 + di 


„ 1 dA 1 dd 2 (71-72) 

» + rfi + 


di 


in (-d,-d 2 ), 


{ C k (-d) = c k (-d 2 ) = 0 , 

whence 


sinh ((d + y)R k ) , 7i 2 , 

Ck{i) =A sinh (RiTi) + V + 


d + d 2 A - ~f 2 d 2 
7i + 


dA ^ 2 ( 71 - 72 ) 

!,+ 7' + —*—' 


di di 

The representation of df& 2 (A, 0) as a Fourier multiplier follows now easily. 

The derivative dh& 2 (A,0) Observing that 

3/ l < MA,0)[h] = —-h + dy{d h w 2 (A,0)[h})\ y= _ d2 , 

«2 

the desired representation for c^d^A, 0) follows by using the expression for dh.w 2 (A,0)[h\ 
determined in the second part of this proof. □ 


I 11 the remaining part we establish the Lemmas A.1-A.6 that provide the justification for 
the streamlines pattern, as seen from a reference frame moving with the wave, as shown 
in Figures 1-3. Because the proofs of Lemmas A.1-A.6 use similar arguments, we present 
herein only the proof for Lemma A.l. For this, it is important to note that because there is 
no time dependence in problem (2.1) (or (2.2)), the particle trajectories and the streamlines 
corresponding to the solutions found in Theorems 3.4-3.6, 4.2, 4.3 coincide with the level 
curves of the corresponding stream function. They are parametrized by solutions of the 
system of ordinary differential equations 


j x' = u — c = ijjy, 
1 y' = v= 


(A.7) 


stagnation points of the flows corresponding to equilibria of (A.7). Hence, our task is to 
determine the level curves of the stream function. The direction of motion of the particles 
along the level curves is determined by the sign of u — c or v. 
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Lemma A.l. Assume that 71 > 72 > 0 and let 

€ (c 3 + a (R)) 2 x c 3 + a (nU)) X c 3 p + a mJM) 

he a solution of ( 2 . 2 ) that is determined by a point (A(s), (f,h)(s)) on the bifurcation curve 
found in Theorem 3.f. Provided that s is small enough, the following assertions are true: 
(i) f > 0 and h! > 0 on (0, L/ 2); 

(ii) d x ;tp 2 < 0 in {(x,y) £ il(f,h) : x £ (0, L/2)} and d y ip 2 > 0 in U(f,h); 

(in) d x ipi < 0 in {(x,y) £ Q(f) : x £ (0,L/2)}; 

(iv) There is a smooth curve {(x,£(x)) : x € [0, L/2]} with —d < £(x) < —d 2 + f(x) for 
all x £ [0, L/2] and satisfying additionally: 

(a) Given x € [0,L/2], it holds that: d y ipi(x,£(x)) = 0, d y ipi(x,y) < 0 for all 
y £ ]-d,£(x)), and d y ipi(x,y) > 0 for all y £ (£(x),-d 2 + f(x)]; 

(b ) £ is strictly decreasing on [0, L/2]; 

(c) The function [x 1-7 ipi(x, ^(x))( is strictly decreasing on [0, L/2]. 


Proof. Since A(0) = Ai £ ( 72 ^ 2 , 71 ^ 1 + 72 ^ 2 ), for small s it holds A(s) £ ( 72 ^ 2 , 71 ^ 1 + 72 ^ 2 )- 
Recalling that 


2 t r 


f(s) = sm\ (Ai)cos ( —x ) + 0(s 2 ), h(s ) = — sm{ (Ai)cos 


2vr 

T : 


+ 0(s 2 ), (A.8) 


with m 22 (Ai) < 0 and m 21 (Ai) > 0, the claim (i) follows by using the same arguments as 
in the proof of [35, Lemma 4.2]. 

For (ii ), we see first that dyip® = liV + Ai > 72 y + 72^2 > 0 in Q 2 . Therefore, d y ip 2 > 0 
in f l(f,h) provided that s is small. Using now (i) and the fact that if is constant on 
<9f l(f, h ) and even with respect to x, it is easy to see that d x tp 2 < 0 on the boundary 
of the set {(x,y) € f l(f,h) : x € ( 0 ,L/ 2 )}. Observing that d x ip 2 (x,h(x)) < 0 for all 
x £ ( 0 , L/ 2 ) and Aip x = 0 in Q(f,h), elliptic maximum principles ensure that d x ip 2 < 0 in 
{(x,y) € Tt(f,h) : x € ( 0 ,L/ 2 )}. The claim (in) is obtained in a similar manner. 

For (iv), we remark that 


dyi’ll = _ d < 0, 


d y^\v=-d? > °> d yy^i > 0 in n i- 


(A.9) 


1 y= — d ’ y + 1 \y=-d2 

Therewith, for small s the function ^\ satisfies the similar inequalities 

dyi’i\ y= _d < 0 , dytpi\ y= _ d2+f > 0 , dyyi/f i > 0 in Q(f). 

Hence, for each x £ [0, L/2], there exists a unique £(x) £ (—d,—d 2 + f(x)) such that 
d y ip\(x, f(x)) = 0. Due to the third inequality in (A.9) we conclude from the implicit 
function theorem that £ is smooth and 

d xy ip i(x,£(x)) + £'(x)d yy ipi(x,£(x)) = 0 for all x € [0,L/ 2 ]. (A.10) 

We are going to determine now the sign of d xy ip \. To this end note that ip\ = uqo ^/ 1 where 

w± £ c 3 +?r(n 1 ) is the unique solution of the problem (2.9), that is w\ := uii(A(s), (/, h)(s)). 
By the chain rule we get 

dif Q _ _ , di Q d\fy 


dxy'lfr 1 — 


77T7F d,m ° +' + dtTJ a * yW1 ° *7' ~ WTJ? BvvW ' ° + 1 


ddi //' ddif . 

+ 7TT7P+’" 1 ° + “ WT7f yyWl ° ’ ■ 

On the other hand we have the following expansion 

wi(A(s),(f,h)s)) = uq(Ai,0) + d\wi(Ai, 0)[A(s) — A x ] + dfw^Ai, 0)[/(s)] + 0(s 2 ) 
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in Observing that 

dyWi o 4 *- 1 = d y i >? + O(s), ) 

dxyW 1 O &J 1 = d xy (d f wi(Ai,0)[f]) + 0(s 2 ), > in C e 2 +" (Ll(f )), 

o <F / 1 = 71 + O(s), J 

a lengthy calculation leads us to 

d xy ^i = — smf (Ai)LiAi sin ( L i+) + °( s2 ) in )■ 

A similar argument to the one used in (i) shows that d xy ipi > 0 in Ll(f) if s > 0 is 
sufficiently small. The latter property together with (A.9) and (A. 10) implies that £' < 0 
in x G (0, L/2). This proves the claim in (b). Since (c) is an obvious consequence of (Hi) 
we have completed the proof. □ 

It follows now readily from Theorem 3.4 and Lemma A.l that the streamline pattern in 
the moving frame for the non laminar solutions found in Theorem 3.4 for 71 > 72 is as in 
Figure 1 (left picture). The next lemma justifies the right picture of Figure 1. 

Lemma A.2. Assume that 71 < 72,72 > 0 and let 

(C€ (qJ+?(R )) 2 x c 3 + a (nU)) x cl+ a (WM) 

be a solution of (2.2) that is determined by a point (A(s), ( f,h)(s )) on the bifurcation curve 
found in Theorem 3-4 ■ Provided that s is small enough, the following assertions are true: 
(i) f > 0 and h' < 0 on (0, L/2) ; 

(ii) dxipi > 0 in {(x,y) G f2(/) : x G (0, L/2)} and d y ipi < 0 in fl(/); 

(Hi) d x if 2 > 0 in in {(x,y) G Q(f,h) : x G (0,L/2)}; 

(iv) There is a smooth curve {(x,£(a:)) : x G [0, L/2]} with — d 2 + f(x) < £(x) < h(x) 
for all x G [0, L/2] and satisfying additionally: 

(a) Given x G [0, L/2], it holds that: d y i/j 2 (x,^(x)) = 0 7 d y f) 2 (x,y) < 0 for all 
y G [~d 2 + f(x),f(x)), and d y if 2 (x,y) > 0 for all y G (£(x),h(x)]] 

(b) £ is strictly decreasing on [0, L/2]; 

(c) The function [x 1 —> if 2 (x,f t (x))\ is strictly increasing on [0, L/2]. 

The next lemma provides a justification for the left picture of Figure 2. 

Lemma A.3. Assume that 72 > 0 , 71 ^ 1 + 72^2 >0 and let 

((f,h),i>iM € (c p 3 + a (M )) 2 X c 3 + a (nU)) x c 3 + a ( ncOT) 

be a solution of (2.2) that is determined by a point (A(s), ( f,h)(s )) on the bifurcation curve 
found in Theorem 3.5. Provided that s is small enough, the following assertions are true: 
(i) f > 0 and h! < 0 on (0, L/2); 

(ii) d x ij) 1 > 0 in {(x,y) G Ll(f) : x G (0, L/2)} and d y i/\ < 0 in fi(/); 

(Hi) d x if 2 > 0 in in {(x,y) G Q(f,h) : x € (0, L/2)}; 

(iv) There is a smooth curve {(a;,£( 2 ;)) : x G [0,L/2]} with —d 2 + f(x) < £(x) < h(x) 
for all x G [0, L/2] and satisfying additionally: 

(a) Given x G [0, L/2], it holds that: d y if 2 (x, f(x)) = 0, d y if 2 (x,y) < 0 for all 
y G [~d 2 + f(x),f(x)), and d y i/j 2 (x,y) > 0 for all y G (£(x),h(x)]; 

(b) £ is strictly decreasing on [0, L/2]; 

(c) The function (x >—>■ if 2 (x,f(x))\ is strictly increasing on [0, L/2]. 

We provide now a justification for the right picture of Figure 2. 
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Lemma A.4. Assume that 72 > 0, 71^1 + 72^2 < 0 and let 

iM € (c p 3 +“(m )) 2 X c 3 + a (nU)) x c 3 p + a mTh)) 

he a solution of (2.2) that is determined by a point (A(s), (f,h)(s)) on the bifurcation curve 
found in Theorem 3.5. Provided that s is small enough, the following assertions are true: 
(i) f > 0 and h' < 0 on (0, L/ 2) ; 

(**) d x if) 1 > 0 in {(x,y) G Q(f) : x G (0,L/2)}; 

(in) d x ip 2 > 0 in in {(x,y) G Tl(f,h) : x G (0, L/2)}; 

(■ iv) There is a smooth curve {(x,£i(x)) : x G [0,L/2]} with —d < fi(x) < —d 2 + f(x) 
for all x G [0, L/2] and satisfying additionally: 

(a) Given x G [0,L/2 ] 7 it holds that: d y ij}\(x,^ i(x)) = 0 7 d y ij}\(x,y) > 0 for all 
y G [~d,f i(x)), and d y ifi(x,y) < 0 for all y G (t >1 (x),-d 2 + f(x)]\ 

(b) £1 is strictly increasing on [0, L/2]; 

(c) The function [x 1-7 i(x,£i(x))] is strictly increasing on [0, L/2]. 

(u) There is a smooth curve {(x,^ 2 (®)) : x G [0,L/2]} with —d 2 + f(x) < &(x) < h(x) 
for all x G [0, L/2] and satisfying additionally: 

(a) Given x G [0, L/2 ] 7 it holds that: d y 'ip 2 (x,£, 2 (x)) = 0, d y ip 2 (x,y) < 0 for all 
y G [~d 2 + f(x),£ 2 (a0), and d y ij> 2 (x,y) > 0 for all y G (&(x),h(x)]; 

(b ) £2 is strictly decreasing on [0, L/2]; 

(c) The function [x i— > if 2 (x, £ 2 ( 2 ))] is strictly increasing on [0, L/2]. 

We consider now the non laminar flows corresponding to the bifurcation solutions found 
in Theorem 3.6 and prove the following result which justifies the left picture of Figure 3. 

Lemma A. 5. Assume that 72 > 0 , 7 idi + 72^2 < 0 and let 

((f,h),if>iM G (qL + r“W ) 2 x c 3 + a (nU)) x c 3 + a (W,jfj) 

be a solution of (2.2) that is determined by a point (A(s), ( f,h)(s )) on the bifurcation curve 
found in Theorem 3.6. Provided that s is small enough, the following assertions are true: 

(■ i) f < 0 and h! < 0 on (0, L/2); 

(ii) d x f )2 < 0 in {(x,y) G Q(f,h) : x G (0, L/2)} and d y i\) 2 < 0 in Q(f,h); 

(in) d x fji < 0 in {(a;,y) G Sl(f) : x G (0,L/2)}; 

(iv) There is a smooth curve {(x,£(x)) : x G [0, L/2]} with —d < £(x) < —d 2 + f(x) for 
all x G [0, L/2] and satisfying additionally: 

(a) Given x G [0, L/2], it holds that: d y ipi(x,£(x)) = 0, d y ifi(x,y) > 0 for all 
y G [~d,£(x)), and dyifi(x,y) < 0 for all y G (£(x),-d 2 + f(x)]\ 

(b) £ is strictly increasing on [0, L/2] / 

(c) The function [x 1 — > i/ji(x,£(x))] is strictly decreasing on [0, L/2]. 

Finally, we have the following result which justifies the right picture of Figure 3. 

Lemma A. 6 . Assume that 72 = 0, 71 < 0, and let 

((f,h),i>iM e (cj£ a (R )) 2 x c 3 + a (nU)) x c 3 + a m7M) 

be a solution of (2.2) that is determined by a point (A(s), (/, h)(s)) on one of the bifurcation 
curves found in Theorems f.2-4-3. Then, the assertions from Lemma A.5 hold verbatim. 
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